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We investigate spectral properties of a 1-dimensional piecewise linear intermittent map, which 
has not only a marginal fixed point but also a singular structure suppressing injections of the orbits 
into neighborhoods of the marginal fixed point. We explicitly derive generalized eigenvalues and 
eigenfunctions of the Frobenius-Perron operator of the map for classes of observables and piecewise 
constant initial densities, and it is found that the Frobenius-Perron operator has two simple real 
eigenvalues 1 and Ad £ ( — 1, 0), and a continuous spectrum on the real line [0, 1]. From these spectral 
properties, we also found that this system exhibits power law decay of correlations. This analytical 
result is found to be in a good agreement with numerical simulations. Moreover, the system can 
be extended to an area-preserving invertible map defined on the unit square. This extended system 
is similar to the baker transformation, but does not satisfy hyperbolicity. A relation between this 
area-preserving map and a billiard system is also discussed. 

PACS numbers: 05.45.-a 



I. INTRODUCTION 

In the past decades, a lot of studies have been de- 
voted to investigations of the relations between micro- 
scopic chaos and non-equilibrium behaviors such as re- 
laxation and transport, and it has been found that mi- 
croscopic chaos plays essential roles in non-equilibrium 
processes [1, 2]. For example, it is well known that, for 
fully chaotic (hyperbolic) systems, correlation functions 
decay exponentially and their decay rates are character- 
ized by the discrete eigenvalues of its Frobenius-Perron 
(FP) operator [3-7]. As one of the examples of the hy- 
perbolic systems that permits detailed calculations, the 
baker transformation has been studied extensively and 
its spectral properties of FP operator are fully under- 
stood [8] . In addition, the baker map is considered as an 
abstract model of chaotic Hamiltonian systems, because 
it has the area-preserving property, which is an universal 
feature of the Poincare map of the Hamiltonian systems 
of 2 degrees of freedom [9]. In fact, similarities of the 
baker map to the Lorentz gas with finite horizon have 
been pointed out [10]. 

In contrast to hyperbolic systems, dynamics in generic 
Hamiltonian systems is more complicated and diverse. 
When a phase space of a Hamiltonian system con- 
sists of integrable (torus) and non-integrable components 
(chaos), power law decay of correlations is frequently ob- 
served [11-16]. Although such kind of systems, i.e., sys- 
tems with mixed type phase spaces, are more generic 
than the integrable or the fully chaotic systems, the theo- 
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retical understanding of their statistical properties is not 
enough; for example the ergodic and mixing properties 
of chaotic components of generic systems are still unclear 
from the theoretical point of view. 

For understanding sub-exponential decay of correla- 
tion functions in dynamical systems, non-hyperbolic 1- 
dimensional maps have been studied by several authors 
[17-26] and they have found power law decay of correla- 
tions in their models. Therefore, it is natural to imag- 
ine connections of these non-hyperbolic maps and mixed 
type Hamiltonian systems, however, extensions of these 
maps to 2-dimensional area-preserving systems are un- 
known. Thus in this paper, we introduce a modified 
version of the 1-dimensional intermittent map studied 
in Ref. [19] and extend it to an area-preserving system. 
This area-preserving map, which is similar to the baker 
transformation, may be considered as an abstract model 
of mixed type Hamiltonian systems. 

Our theoretical treatment is mainly based on Ref. [19], 
where a piecewise linear version of the Pomeau- 
Manneville map [27, 28] is proposed and its general- 
ized spectral properties of the FP operator in a sense 
of Refs. [29, 30] have been elucidated. Their model has 
a marginal fixed point and exhibits power law decay of 
correlation which they have found to be an outcome of a 
continuous spectrum of the FP operator. In addition to 
a marginal fixed point, the piecewise linear map studied 
in the present paper has a singular structure, which sup- 
presses injections of the orbits into neighborhoods of the 
marginal fixed point. Due to this property, the uniform 
density is invariant under time evolution and the map 
can be extended to an area-preserving map on the unit 
square. And it is shown that generalized eigenvalues of 
the FP operator consists of two simple real eigenvalues 1 
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FIG. 1: The piecewise linear map <f){x) (solid line) for b = 0.5 
and f3 = 1.3. The circles indicate endpoints of the straight 
line segments. 



and Xd S (—1,0), and a continuous spectrum on the real 
interval [0, 1]. It is also shown that correlation functions 
exhibit power law decay due to the continuous spectrum. 

This paper is organized as follows. In Sec. II, we in- 
troduce the piecewise linear map, and define the FP op- 
erator and functional spaces of observables and initial 
densities. In Sec. Ill, we derive the spectral decomposi- 
tion of the FP operator. In Sec. IV, long time behaviors 
are analyzed and some numerical results are displayed. 
We also discuss the extension of our model to an area- 
preserving invertible map in Sec. IV. Section V is devoted 
to summary and remarks that include comments about 
differences between our model and the one in Ref. [19], 
and about similarities to a billiard system. 



II. PIECEWISE LINEAR MAP 

A. Definition of the map 

The dynamical system we consider in this paper is the 
piecewise linear map shown in Fig. 1. This map (t>{x) : 
[0, 1] [0, 1] consists of two parts: the left < x < 6 
and the right part 6 < .x < 1. Each part is formed by 
an infinite number of straight line segments; the circles 
plotted in Fig. 1 indicate endpoints of these segments. 
The map ^(a;) is defined by 

forage [4-, (fc = l,2,...), 
forxe (fc=l,2,...). 



(i){x) = < 



(1) 



In this definition, ^j, (fc = 0, 1, • • •) represent the horizon- 
tal coordinates of the endpoints of the segments on the 



left part {x < b) and they are defined as 
Co" = b, 



^k-l ~ 



for fc = 1, 2, • • • . 



(2) 



where /3 > 1 is a parameter and = J2^=i ^/n^ is 
the Riemann zeta function. And r]'^ is the slope of the 
fc-th segment and defined as 



i 



Ck 



k-2 Sfe-l 

1-6 



C(/3) for k = l, 



k-1 



(3) 



(4) 



for fc = 2, 3, • • • , 



where wc; dcifinc — 1 for convenience. This definition 
of the left part (0 < a; < 6) is the same as that of the 
piecewise linear Pomeau-Manneville map proposed in the 
Ref. [19]. The map ({){x) can be approximated as (t>{x) ~ 
X + Cx^/^l^'^^ when x 0+, where C is a constant. 
Thus the origin a; = is a marginal fixed point. 

In the same way, (fc = 0, 1, • • •) are the horizontal 
coordinates of the endpoints of the segments on the right 
part and defined as 



^+ = 6 1 + 



1 



C(/3) 
b 

for fc = 2, 3, • • • 



^0+ = 1> 



e+ - f + - 



(fc - ly fc/3 



(5) 



and rf^ is the slope of the fc-th segment and defined as 



4 



= < 



1 



l-6[l + l/C(/3)] 
fc/5 



I fc/3-(fc-l) 



for fc = 1, 



for fc = 2,3,-. 



(6) 



(7) 



This is the definition of the right part of the map (t){x), 
b < X < I. This part is different from the Pomeau- 
Manneville type maps, and is similar to a map proposed 
by Artuso and Cristadoro [18]. 4'{x) behaves as (j){x) ~ 
{x — b)^^~^^/^ when x b+. Therefore the derivative 
4>'{x) of the map is divergent at a; = 
We assume that ryj*" > 0, i.e.. 



b< 



m 



(8) 



Note that the uniform density on the interval [0, 1] is 
invariant under the time evolution of this map because 
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the relation l/?7^ + ^/Vk = 1 is satisfied for fc = 1, 2, ■ • •. 
Figure 1 shows the shape of the map 4>{x) for (3 — 1.3 and 
h = 0.5. There is a singular structure near x = b, which 
suppresses injections of the orbits into neighborhoods of 
the marginal fixed point x = 0. This system can be easily 
extended to a 2-dimensional area-preserving map, whose 
dynamics of the expanding direction is given by the map 
(pix). This will be discussed in Sec. IV. 



B. FP operator and functional spaces 

The FP operator P and its adjoint P* are defined by 



Pp{x) = [ dyS{x-(f){y))p{y), 
Jo 



P*A{x) = Ai<j>{x)), 



(9) 



(10) 



respectively. We also define an inner product {A, p) as the 
average of an observable A{x) with respect to a density 
P{x), 



{A,p) = / dxA(x)p{x). 
Jo 



(11) 



Then, the average of A{x) at time t with respect to an 
initial density p{x) is given by {A,P*p) = {P*^A,p). 

Let us consider an observable A{x) such as the inequal- 
ity 



\A{x) - ao - aix\ < Kx^^/^^-^^ 



(12) 



holds for some positive constant K, where the constants 
ao and ai satisfy ao = ^(0) and ai = ^'(0), respectively 
[19]. This function is bounded on [0, 1] and smooth near 
the origin. And we define a set Xq as the functional 
space which consists of such observables. This functional 
space is invariant under the action of the adjoint of the 
FP operator P*, namely, P*A{x) G Xq if A{x) e Xq. 
If we define the norm as 



P(a;)||o = |ao| + |ai|+sup|A(a;)| 

X 

\A{:r) - «() - r/i,c| 



+ sup- 



xl^/i0-^) 



(13) 



then this functional space becomes a Banach space with 
respect to this norm. Note that this space is dense in the 
Hilbert space ^^[0, 1] of the square integrable functions 
on [0,1]. 

Furthermore we restrict initial densities to be piecewise 
constant [24], 

p{x)=pk if [4-,4-_i) (fc = 0,l,2,-..). (14) 
We also assume the following properties for k = 1, 2, • ■ •, 



k + l 



with ^|p,|^'<+oo, (15) 



where > 1 is a constant. We also assume the normal- 
ization condition 



ES + (1-^)^o = 1. 



(16) 



In the following sections, we use this condition [Eq. (16)] 
only for clarity of exposition. But it is not essential and 
almost the same result can be obtained without it. 

We define a set Xn as the functional space which 
consists of such densities. This hmctional space is in- 
variant under the action of the FP operator; namely, 
Pp G Xd if p € Xd- And if the norm of this space is 
defined as 



(17) 



1=0 



this functional space becomes a Banach space. The func- 
tional space Xu is not dense in L^[0, 1]. 

In the above definition for initial densities, we assume 
that the initial densities are constant on the interval [6, 1), 
i.e., the right part of the map. Although it seems to be 
a strong restriction and it is possible to extend to the 
densities piecewise constant also in the right part, this 
extension does not make any changes for the long time 
behaviors because relaxation to a constant density in the 
right part finishes by only one iteration of the map (l>{x). 



III. SPECTRAL ANALYSIS 

The purpose of this section is to derive a spectral de- 
composition of the average {A,P*p). First, we derive 
the matrix elements of the resolvent operator of P, then 
the average {A, P^p) is obtained by an integral transfor- 
mation of the matrix elements of the resolvent. Finally, 
deforming the integration path, we have the spectral de- 
composition. 



A. Matrix elements of the resolvent operator 

Let us define the matrix elements of the resolvent of 
the FP operator P as, 



(-.;^'')^§j^(-.^V) 



= ^[pkB,;{z)+poB+{z)\, (18) 
fe=i 

where B^{z) are defined below [see Eq. (21)]. Let us 
rewrite {A, P*p) as 

(a,PV) = ^ dxAo^*{x)p{x) 

OO 

= EMfe(*)+^o^fe(*)]' (19) 



fe=i 
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where B'^{t) are defined for = 1, 2, • • • by 



B. 



dxA o (j)\x). 



Then Bj, (z) are defined for fc = 1, 2, • • ■ by 

- Bl' 



6±/ N _ B^jt) 



(21) 



From Eq. (20), we have the following recursion rela- 
tions for B'^{t): 



^ oc 

Bt{t+1) = —Y,B+{t), 
^1 k=i 



+ = Z±Bk-iit) for fc = 2,3,.--. (23) 



(20) 

From Eqs. (21)-(23), we have the recursion relations for 



(24) 



— B^_i(z) for fc = 2,3,---. 
^^fe (25) 



Using these relations recursively, the following equation 
(22) can be derived for fc = 1, 2, ■ • •: 

I 



m=0 



k 



zkk0 (1 - 6)C(/3) 



(26) 



Furthermore, from Eqs. (25) and (26), we obtain for A: = 2, 3, 



Bim , 1 



k-2 



i(0) fk-m-1 



.m+2 



.m=0 



k-1 

r 



z^{k - ly (1 - 6)C(/3) 



m=l 



(27) 



Summing up Eq. (24) and Eqs. (27), we have 

fe=i ^ 
where Z(^;) and '^{z) are defined by 



and 



Z{z) 



r?+ (1 - 6)C(/3) ^ ry^+i^^fc/^ 



E;;t- 



$(.) ^ E^fe(o)+EE 



BiiO) 



fe=i 



^ ^ 77+ , , , \ni + I 

1 rn=0 \ 



(28) 
(29) 

r- 

(30) 



H(z) = B^m{z)+j2Biiorz'-' 



*(2) - Po 



1=2 



Pk 

(1 - ^')C(/3) t{ '"f"^ 



■E 



(33) 



respectively. These functions ^'(z) and ^(z) are also ab- 
solutely convergent and analytic for \z\ > 1. 



B. Analytic properties of individual functions 

For r > 1, we obtain the average {A, pt) by an integral 



The functions Z{z) and <I>(z) arc absolutely convergent 

for \z\ > 1 and thus analytic there. From Eqs. (26), and , ^ "^IT"" ""^ "'^^7^'" ^'^"oi^: ""'r^n ""'^^ 

/r,o\ • 1 r.Li. i /i oM transformation of the resolvent Eq. (31) as foUows: 

(28)-(30), the matrix elements of the resolvent [Eq. (18)] l -i \ /j 

can be rewritten as 



A, 



z-P' 



^(z)$(z) 



EE 

;=i fe=; 



^k-l+l 



*'-'«^''<'-*'««[H(.)-Sr(0)«l,(31) 



|z|=r- 



dz f 
2m 



vl>(z)$(z) , (1^&)CM5(^)' 



Z(z) 



2-(z) 6 
where we define the functions 4'(z) and S(z) as 

h ^ Pk 



(34) 



'J'(2:) =Po + 



(1 - 6)C(/3) z'^ 



(32) 



where the integration path is taken in the counter clock- 
wise direction. In order to deform this integration path 

into the unit disk 1^1 < 1 and derive the spectral decom- 
position, we study the analytic properties of the functions 
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in the integrand of Eq. (34) in this subsection. With the we have analytic continuations of the functions Z{z), 
help of the identity <&(z) and \E'(-2) into the unit disk \z\ < 1. For example, 



1 1 



k0 r(/3) 



POO 

/ ds s^-^e-^', (35) 

^0 



^{z) can be analytically continued as follows: 



dssP-^e-^\l - 6-^)^(2-^6-^)™ 



m=0 



k=l ^' ' 1 = 1 •'^ 



(36) 



(37) 



where these calculations are justified because the convergence of the summation is uniform in s for \z\ > 1. Note that 
each integral in the right-hand side (rhs) of Eq. (37) is analytic except for the real interval [0, 1]. Then the function 
$(z) expressed as Eq. (37) is analytic except for the cut [0, 1], because the second infinite sum in Eq. (37) is absolutely 
convergent, i.e.. 



1 



B-{or 



ds 



z — e~ 



< 



bsnp\A{x)\ 

X 



(38) 



where d{z, [0, 1]) is the distance between z and the real interval [0, 1]. 

In the same way, we obtain analytical continuations of the functions Z{z) and "^{z): 



Z{z) = {z-\) 
^(z) = po 



(1 - 6)C(/3)r(/?) Jo 



ds 



z — e 



E 



pi 



oo ^/3-lg-(i+l)s 

ds . 

z — e ^ 



(39) 



(40) 



These expressions for the functions Z{z) and ^{z) are also analytic except for the cut [0, 1]. And the function '^{z) 
can be also analytically continued into the unit disk \z\ < 1 in the same way. Here, however, we rewrite the function 
2(2) in the following form: 



E{z) = ^{z)J2l^Bf{0)z'-' -J2l^Bi{0)z'-^ 



^ V 

i^-b)mti 



Pk 



(41) 



The infinite sums in this expression of the function '^{z) are absolutely convergent for \z\ < 1. 



r 



Let us consider the zeros of the function Z{z). First 
wo define a function Q,{z) as Q.{z) = Z{z)/{z — 1). And 
if Im{z) + 0, then 



Im Q.{z) 



6Im(z) 



(1 - 6)C(/3)r(/3) Jo k - e 



ds 



s^-^e-' 



^0. 



(42) 

And also if Imz = and Rez > 1, then f7(z) > 0, because 
z — 6~* > 0. Thus there is not zero of Z{z) in these 
regions. 

Next if Imz = and Rez < 0, we have 



n'{z) = 



(1 - b)mm Jo - e-«) 



ds- 



< 0, (43) 



and 



> 1 



= 1 



{i-b)mm Jo 

b 

{i-b)mm 

b 



f 

Jo 

f 

Jo 



ds 



s^-^ 
+ 1 

O0-1 



ds- 



(1 - bKW 



> 



(44) 



The last inequality holds because of the Eq. (8). And it 
is easy to see that Cl{z) — > —00 as ^; — > 0. Therefore, on 
the real interval [—1,0], fl{z) has the unique zero of 
order 1. 
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C. Spectral decomposition 

Prom the results of the last subsection, the contour of 
the integration of Eq. (34) can be deformed into the unit 
disk \z\ < 1 as 



Z{z) 



+ lim 



dz 



*(z)<I>(z) ^ (l-6)C(/3) ^^,^- 



Z{z) ' b 
where the integration path C is defined by 

C = {z\\z\=v, \z-l\=v 

or z = X±iO {r] < X <1 - T])} . 



,(45) 



(46) 



The first term of the rhs of Eq. (45) corresponds to the 
simple pole at ^; = Ad and thus can be calculated by the 
residue theorem and Eq. (45) is rewritten as 



(Ad - i)n'(Ad) 



+ lim 



■^{z)^{z) {i-b)m 

Z{z) b 



E{z) 



.(47) 



Let us consider the integral of rhs of Eq. (47). Because 
of the facts (see the appendix) 



Jo 



dxA{x), 



(51) 



where we have used the normalization condition Eq. (16). 
The rhs of Eq. (52) is the average of A{x) with respect 
to the invariant density which is uniform for the map 
(pix); therefore Eq. (52) gives the average value of the 
invariant state. 

Finally wc have to evaluate the integral along the cut. 
For this purpose, we define the new functions Cl{X), fx{x) 
and ^[{x) for < A < 1 as follows, 



f2(A) = 1 + 



dsV- 



l-b 



(1 - 6)C(/?)r(/3) Jo A - e 



^ oo 

xo(a;) + |^5^?''xi(a;) 



(52) 



1=1 



L 



X / dsV- 



X-e- 



°° g/3-lg-(i+l)s 



A-e- 



(53) 



(54) 



where V means the Cauchy's principle value and xii?^) is 
defined, for Z = 0, 1, ■ • •, by 



and 



lim z°''E.{z) = 0, 

2— >0 



lim ^" ^ J,, } ' = 0, 



Z{z) 



(48) 



(49) 



for > 0, the contribution from the integration around 

the origin, {|z| = rj}, vanishes. 

Similarly, we obtain (sec the appendix) 



\im{z-l)E{z)=0, 

Z—>1 



(50) 



and then the contribution from the integral around z = 
1, {|z — 1| = T]}, vanishes for the second term of the 
integrand. On the other hand, we have 



hm / 

''^1 J|0-l|=77 



dzz 



,*(£)£(£) 

Z{z) 



^, Hz)^{z) 
lim (2; — 1 — , 
^^1^ ' Z{z) 



1 for < a; < ^i-i, 



otherwise. 



(55) 



Note that wc have defined = 1 in Sec. II. These 
functions are related to the real parts of the functions 
{Z{z),^{z),^{z)} near the cut: ReZ(A ± iO) = (A - 
l)f2(A), Rc$(A ± iO) = {I - b) JldxA{x)fx{x), and 
Re«'(A ± iQ) = po + &/(l - b) YZo pM^)- 

It can also be shown for the imaginary parts of the 
functions {Z{z),^{z),'^{z)} near the cut that 



Im$(A - ?:0) _ (1 - b)C{f3) 
\mZ{X + «0) ~ b 



Y,l^B-{^)^-\ (56) 



Im*(A - iO) 
ImZ(A + iO) 



00 



(57) 



/' 

JO 



dxA{x) 



(1 - h)po + 



E 



Pk 



Using these functions, we can calculate the integral 
along the cut as follows [19]: 
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lim 



dz 



"J^O Jc\{\z\=v U \z-l\=v} ^""^ 



^ rfA f Im [Z(A + zO)$(A - iO)] Im [Z{X + 20)1'(A - zO)] 
\Z{X + iO)\nm Z{X + iO) 



(1 - b)amp) Jo 



1 Im [Z{X + iO)$(A - lO)] A* (1 - A) (log i 
uA- 



ImZ(A + iO) |Z(A + iO)|2 

= [\x{A,Fx)X'{Fx,p), 
Jo 

where we define the hnear functional {A, Fx) of the observables A{x) as 



Im [Z(A + iO)*(A-iO)] 
Im Z{X + iO) 



{A, Fx) = N{X) f dx{A{x) - A(0)} 
Jo 

The function A''(A) in Eq. (61) is given by 

JV(A) = 



(58) 
(59) 
(60) 

(61) 



C(/3)r(/3)(i-b) 



(l-A)<^f22(A) + 



(log 



(62) 



We also define in Eq. (60) the linear functional {Fx , p) of 
the initial densities p{x) as 



(68) 



(/'A,p) = (l-6)po + ^>^Pi 



z=o 



^i(A) 



1-6 



f7(A)A' 



(63) 

Similarly, the linear functional associated to the eigen- 
values z — Xd and z = 1, that have been obtained in 
Eqs. (47) and (52), can be expressed as 



(A^k) 
{Fd,p) 



1 



{Xd - i)i^'(Ad) y„ 

oo 

(1 - b)po + h^piVi{Xd) 



dxA{x)fxAx), (64) 
(65) 



and 



{A Fin) 

{Fin,p) 



= / dxA{ 
Jo 



x), 



(66) 
(67) 



respectively. It is easy to check that {A, Fx) = {A, Fd) = 
if A{x) is a constant function and that {Fx, p) = 
{Fd,p)=Oif p{x) = l. 

Using these linear functionals, we can spectrally de- 
compose the average {A,pt) [Eq. (47)] at time t = 
0, 1, 2, • • • as, 

{A,P'p) = {A,Fin){Fin,p) + Xd'{A,Fd){Fd,p) 



C(/3)r(/3)(l-6) 

I 

+ [' dXX'{A,Fx){Fx,p). 
Jo 

The left eigenfunctions {.Fin, -Fd, -Fa} satisfy the relations 

{A,PFin) = (P*A-Fln) = {A, Fin), 
{A,PFd) = {P*A,Fd) = Xd{A,Fd), (69) 
{A,PFx) = {P*A,Fx) = X{A,Fx). 

Therefore {Fin, -Fa j -Fa} are eigenfunctions of the FP op- 
erator in a generalized sense [19, 24, 29, 30]. On the 
other hand, the right eigenfunctions {Fin, -Fd, -Fa} satisfy 
the relations 

(F*Fin,p) = (Fi„,Pp) = (Fi„,p), 
(F*Fd,p) = (Fd,Fp) = Ad(Fd,p), (70) 
{P*Fx,p) = {Fx,Pp) = X{Fx,p). 

Therefore {Fin,Fd,-FA} are eigenfunctions of the adjoint 
of the FP operator in a generalized sense. 

Thus wc obtain the spectral decomposition Eq. (68), 
where the spectrum consists of two discrete eigenvalues 
1 and Ad, and a continuous spectrum on [0,1]. And 
Eq. (68) for t = shows that these eigenfunctions are 
complete. 



8 




FIG. 2: (a) Power spectral densities S{f) of time series A{x{t)) (in log-log form) for five different values of the system parameter 
13 : (3 =1.2 (the solid line), 1.4 (dotted), 1.6 (broken), 1.8 (dashed), 2.0 (dashed-dotted) , 2.2 (dashed- two-dotted) . These PSDs 
are obtained by averaging over 20000 initial conditions uniformly distributed in [0, 1]. (b) The scaling exponent 7 of the PSD 
5'(/) ~ 1//'' as a function of j3. The circles are the numerical results obtained by least square fitting in the low frequency 
region (below / = 10"'') of the PSDs S{f)\ and the dashed line is the theoretical prediction [Eq. (73)]. 



IV. LONG TIME BEHAVIORS AND AN 
AREA-PRESERVING EXTENSION 

A. Long time behaviors 



where K' is a constant. Eq. (72) shows that the corre- 
lation functions decay algebraically. From Eq. (72), it 
is found that the power spectral density (PSD) S{f) be- 
haves as [25] 



In the previous section, we derive the spectral decom- 
position of the average [A, P*p) and found that there is 
a continuous spectrum. Since long time behaviors are 
controlled by eigenvalues whose absolute values are close 
to 1, we consider the limit A ^ 1 for the continuous 
spectrum. The cigcnstatc associated with the eigenvalue 
Xd does not contribute to long time behaviors, because 
this state decays exponentially. Under the assumption 
Eq. (12), the leading term of the left eigenstate {A, Fx) 
when A ~ 1 is given by 

{A, Fx) K{1 - A)'3-2 dx[A{x) - A{0)] (71) 
Jo 

where if is a constant. Similarly, we obtain the leading 
term for the right eigenstate {Fx, p), 

00 CXD 

{Fx, p) ~ (#in, p) - ^Pl = l-^Pl 



1=0 



1=0 



as A ~ 1. Using these facts and Eq. (68), we have for 
t ^ 00 



{A,P'p) 



dxA{x) 



+ 



dx[A{x) - A{0)], (72) 



S{f) ~ <^ 



f 1 




for 




' |log/| 


for 


const. 


for 



for 1<(3<2, 



(73) 



Figure 2(a) shows PSDs S{f) of time series A{x{t)), 
where the observable A{x) is a step function defined on 
X € [0, 1) as 



A{x) 



-1 for xe [0,1/2) 
1 for xe [1/2,1). 



(74) 



And x{t) is produced by the successive iterations of the 
map 4'{x): x{t -|- 1) = (f){x{t)). Each PSD is obtained 
by averaging over 20000 initial conditions uniformly dis- 
tributed in [0, 1]. 

In Fig. 2(a), the PSDs for /3 < 2 exhibit clear l/p 
scalings in low frequency regions. Figure 2 displays this 
scaling exponent 7 of the PSD S{f) ^ l/.f as a function 
of the system parameter /3. And we show the theoretical 
prediction derived in the above [Eq. (73)] by the dashed 
line in Fig. 2(b). Obviously, the numerical results show 
a good agreement with the theoretical prediction. 
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FIG. 3: The area-preserving map 'tp{x,y) defined by Eq. (75) for h = 0.5 and (5 = 1.3. The domains with labels {A, B, C, a, b, 
c} in the left cell (I) arc mapped into the domains with the same labels in the right (II), respectively. Each domain is uniformly 
stretched in the horizontal direction and uniformly squeezed in the vertical direction. There are infinite number of such domains 
and each one is mapped in a similar way. The partition is displayed only for the region {x G [Ci5!5(7]} ^'^'d {x G [Ce^j^iTl} 
the left cell (I); only for x G Krs,^^] in the right (II). The other regions are not displayed because the structures are too fine 
to see. 



B. Area-preserving extension 

As stated in the introduction, the 1 dimensional map 
investigated in the present paper can be extended to an 

arca-prcscrving 2-dimcnsional transformation defined on 



the unit square; this extension 
defined as 



[0,1)2 ^ [0,1)2 



pieces like 

oo 

[0,1)2 ^ \J{{x,y) I a;e [4-,4-_i),t/ G [0,1)} 

oo 

uU{(x,y) |xe[4+,4+_i),2/G[0,l)}. (76) 



fc=l 



% (a; -4) + 4-1' ^ 



for a; e [4 



\ 'Ik 'Ik 



(75) 



% Vk 
fora;e[4+,4+_i). 



where k = 1,2, Obviously, the transformation for 
the horizontal coordinate x is the same as the map (j){x) 
defined by Eq. (1) and docs not depend on the verti- 
cal coordinate y. This relation between 1-dimensional 
map (f){x) and its area-preserving extension xj}{x,y) is 
the same as that between the Bernoulli and the baker 
transformations. Note that the map ip{x,y) is area- 
preserving, because the Jacobian of this map equals to 
1 everywhere. The phase space (i.e., the unit square) of 
this area-preserving map can be partitioned into infinite 



Sec Fig. 3(1). Each piece of this partition {{x,y) \ x £ 
[Cfc^,CfeLi),y G [0, 1)} is uniformly stretched in the hori- 
zontal direction, and uniformly squeezed in the vertical 
direction. The left pieces {{x,y) | x £ [CfeiCfe-i))^ € 
[0,1)} are mapped to the bottom part of the unit cell 
and the right pieces {{x,y) \ x G [^^ G [0, 1)} to 

the upper (Fig. 3(11)). 

Note that although our theoretical and numerical re- 
sults arc for the 1-dimensional map 4>{x), these results 
are also true for this area-preserving map if an observ- 
able does not depend on the vertical coordinate y, namely 
A{x,y) = A{x), because the horizontal coordinate x of 
this area-preserving map is transformed by <p{x) and in- 
dependent of y. Therefore this area-preserving extension 
of (l){x) has also long time correlations with power law 
decay. This fact contrasts to the results for the baker 
transformation, which exhibits exponential decay of cor- 
relation functions. 

Here let us define some terms for later discussions. We 
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define the fc-th escape domain "Dk as 

-Dk = {{x,y) I X e [4--i.4"-2)>y e [0, l)} (77) 

for k = 1,2, • • •. For the points in fc-th escape domain 
Dfe, it takes (k — 1) times of the mappings ijj" to escape 
from the left part a; < 6 to the right x > b. We also define 
the fc-th injection domain as 

= Wx,y) I X G G [0, 1)} (78) 

for k ~ 1,2,---. The injection domains are the upper 
part of the unit square, which are displayed in Fig. 3(11) 
by the labels {A,B,C,- ■ ■}. We denote the areas, namely 
the Lebesgue measures, of the k-th escape and injection 
domains as and 5'^", respectively. These areas obey 
scaling laws Sk ~ l/fc'^ and 4" ~ l/kf^+\ This latter 
power law gives the escape time distribution [31]. 

V. SUMMARY AND REMARKS 

In this paper, we have introduced a piecewise linear 
map and analyzed its spectral properties. We have de- 
rived the generalized eigenfunctions and eigenvalues ex- 
plicitly for classes of observables and piecewise constant 
initial densities. Our model is a modified version of the 
map analyzed in Ref. [19]. A main difference of these two 
models is the normalizability of invariant densities. The 
invariant density of the model investigated in Ref. [19] is 
not normalizable for a parameter region. This is a typ- 
ical property of dynamical systems with marginal fixed 
points [19, 20, 26], and is caused by divergence of invari- 
ant density at marginal fixed points. 

On the other hand, the uniform density is invariant 
for the map 4>{x) discussed in the present paper; there- 
fore the invariant density is normalizable for any values 
of the system parameters, even though our system has 
also a marginal fixed point. This is because the present 
model has the mechanism suppressing injections of the 
orbits into neighborhoods of the marginal fixed point and 
this property prevents divergences of the invariant den- 
sity at the marginal fixed point. As a consequence of the 
normalizability, the present model does not exhibit non- 
stationarity, which is generically observed in maps with 
marginal fixed points [19, 20, 26]. 

The spectral properties of the present model is simi- 
lar to those of Ref. [19] in the locations of the discrete 
and the continuous spectra. There are two simple eigen- 
values 1 and Xd e (—1, 0); the former corresponds to 
the invariant eigenstate and the latter to the oscillating 
one. The eigenstate associated to A^, however, does not 
contribute to the long time behaviors of the correlation 
functions because it decays exponentially fast. There is 
also the continuous spectrum on the real interval [0, 1]; 
this continuous spectrum leads to power law decay of cor- 
relation functions. We have confirmed a good agreement 
between the theoretical prediction and the numerical re- 
sult for scaling behaviors of the PSD S{f) ~ I//''. 



Furthermore, the piecewise linear map (l){x) has been 
extended to an area-preserving invertible map on the unit 
square. In contrast to the baker transformation, which 
is hyperbolic and shows exponential decay of correla- 
tion functions, our model is non-hyperbolic and displays 
power law decay of correlations. 

As is well known, the mixed type Hamiltonian systems 
often exhibit power law decay of correlation functions. 
The area-preserving map ^p{x, y) introduced in this pa- 
per may be considered as an abstract model of the mixed 
type Hamiltonian systems in the following sense. Insta- 
bilities of the orbits of the map ■!/'(x, y) [Eq.(75)] is weak 
in neighborhoods of the line a; = 0, and the escape time 
from the left part a; < 6 to the right part x > b diverges as 
a; ^ 0. In other words, the orbits stick to the line a; = 
for long times. This property seems to be similar to dy- 
namics of Hamiltonian systems near torus, cantorus, and 
marginally unstable periodic orbits, where chaotic orbits 
stick for long times. 

And, in fact, similar dynamics is observed in a Poincare 
map of the mushroom billiard [31], which has been pro- 
posed recently as a model of mixed type systems with 
sharply divided phase spaces [32-36]. In Ref. [31], it is 
found that an infinite partition can be constructed on a 
Poincare surface using escape times from neighborhoods 
of the outermost tori, and that the area of the escape 
domains and the injection domains obey the scaling re- 
lations Vk ^ 1/k^ and V™ ^ 1/fc^, respectively. These 
relations correspond to the case /3 = 2 of the present 
model. Note that a correlation function of the Poincare 
map of the mushroom billiard exhibits power law decay 
C(n) ~ 1/n [31], and this is consistent with the analyt- 
ical result of the present paper. This relation between 
the map 'ij){x,y) and a billiard system is similar to that 
of the baker map and the Lorentz gas [10]. 

Since the map ^{x,y) is an elementary model of con- 
servative systems and can be treated analytically to some 
extent, this system may be important for understanding 
relationships between non-equilibrium phenomena, such 
as relaxation and transport, and underlying reversible 
dynamics; this is a fundamental problem in dynamical 
system theory and statistical mechanics [1, 2]. 
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APPENDIX: ASYMPTOTIC BEHAVIORS 
1. When \z\ 



In this appendix, we show the inequality, 





I 

Jo 



ds- 



<C(log^ 



(A.l) 



as \z\ and argz G [0, 2tt] is fixed, where C is a positive 
constant. From this inequality, the Eqs. (48) and (49) can 
be derived. Let z = x + iy in this subsection. 

First, let us assume arg^ e [37r/4, 57r/4]. Splitting the 
integral into two pieces, we have 



ds 



< 



/•log ■ 

Jo 



dss^- 



+ 



(A.2) 



Apparently, the first term of the rhs of Eq. (A. 2) has an 
upper bound Ci (— log |2;|)'' for some positive constant 
Ci. On the other hand, the second term has an upper 



bound C( (— log |z|)''~"'^. This is obtained by using an 
asymptotic expansion of the incomplete gamma function 
(see e.g., Rcf. [37]). Thus, in this case the inequality 
Eq. (A.l) holds. 

Second, we consider the case for argz € [7r/4, 37r/4] or 
arg.2 G [57r/4, 77r/4]. Using similar calculations, we have 



L 



ds-. 

\z — e~ 



< llog^ 

\y\ 



ds — , 
^/His) 



+ 



POO 

/ ds 



s^-'e- 



\y\ 



-, (A.3) 



where we define i?(s) asH{s) = (a;/|y|^ - 1)^ + lypfi"''). 
It is easy to check that H{s) > 1/2. Therefor, the 
first term has an upper bound C2 (— log jzl)''. By us- 
ing the asymptotic expansion of the incomplete Gamma 
function, we have an upper bound for the second term: 
C^(-log|z|)'^"\ Thus, this case also satisfies the in- 
equality Eq. (A.l). 

Finally, we consider the case aigz G [0, ■7r/4] or aigz e 
[77r/4, 2n]. We split the integral as 



f 

Jo 



ds 



z — e~ 



'G[x-t) G{x + t) 



t + iy 



t-iy 



] + [ dt 

/ J3x/2 



G{t) 
z-t' 



(A.4) 



r 



where we define G{t) as G{t) = {\og(l/t)Y-K By using 
similar techniques as above, it can be shown that the ab- 
solute values of the first and the third terms of the rhs 
of Eq. (A.4) have an upper bound C3 (— log a;)''. For the 
imaginary part of the second term, it is easy to derive an 
upper bound of its absolute value: C3 (— loga;)'^"^. On 
the other hand, for the real part, we have an upper bound 
C3(— logx)^ ^, which can be derived through integra- 
tion by parts after changing the variables as t' = t/x. 
This completes the proof of Eq. (A.l). 



2. When \z\ 1 

The property Eq. (50) is also derived in the same way. 
We briefly mention about the derivation. Let us begin 
with an analytic continuation of '^{z) 



S(z) = B^{Q)^{z) 

00 00 

3=0 1=2 ^ -^0 



ds 



S/3-lg-O>0s 



z — e~ 



(A.5) 



We consider the second term. The absolute value of the 
second term has an upper bound 



const. X 



00 „c 

i=o -^0 



ds 



g/3-lg-04-2)s 



')(1 



(A.6) 



Therefore, we show that for A: = 1, 2, • • ■, and 1 < /3 < 2, 

nl 



ds-. — 
[z ■ 



<G'\l-zf-^, (A.7) 



as z — > 1 and arg(z— 1) G [0, 2tt) with a constant C . Note 
that the integral in Eq. (A.6) from 1 to 00 is convergent. 
For /3 > 2, it can be analyzed in a similar way but the rhs 
of Eq. (A.7) should be changed to a constant (/? > 2) or 
a log correction — log |1 — 2:| {P = 2). Let z = 1 + x + iy 
in the following. 

First, let us assume arg^; e [7r/4, 37r/4] or aigz € 
[57r/4, 77r/4], then we have 



Jo 



ds- 



< e 



Jo 



0/3-2 



ds- 



-, (A.5 



where we have used 1 — e ^ > s/e for s G [0, 1]. The rhs 
of Eq. (A. 8) can be estimated, by splitting the integral. 
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as 



1 gfi-2 

ds-. 

\z-e- 



< 



2e|y| 



«/3-2 



ds- 



\y\ 



J2e\i 



2e / c^ss^-^ 

'2e|j/| (A.9) 



where we have used 1— e ^ > s/e again. It is obvious that 
the rhs of Eq. (A.9) is less than C^|l - z\'^-'^. Therefore 



the inequality Eq. (A. 7) is satisfied in this case. 

Second, the case for argz G [0, 7r/4] or argz S 
[77r/4, 2-k) can be analyzed in the same way as the first 
case (But, in this case, split the integral in terms of x 
instead of y.). Thus we omit the detail. 

Finally, when arg^; e [37r/4, 57r/4], we have 



J 



Jo 



ds 



-lf,-ks 



{z - e-«)(l - e-*) 



-3e/2 p,^. 

dt- ^ ' 



1/e 



Z-t 



.e/2 

/ dt 
Jo 



F{l-e-t) F{l-e + t) 



t + iy 



t-iy 



r 



dt 



l-e/2 



m 

z-t'- 



(A.IO) 



where we define e > as e = and F{t) as F{t) = 
{log(l/t)}'^-it'=~V(l - t)- The first and the third terms 
of the rhs can be estimated in the same way as Eq. (A.9), 
and we obtain an upper bound for their absolute values as 
Csjl — z\^^'^. For the imaginary part of the second term, 
we have easily a bound of its absolute value: C'^\l—z\^~'^. 
On the other hand, for the real part, we also have a bound 



C^'ll 



^ through integration by parts after changing 
the variables ast' = t/e. Thus we complete the proof of 
Eq. (A.7). 

From Eq. (A.7), it can be shown that ^{z) \1>(1) < 
00, as z ^ 1. Thus the first term of Eq. (A. 5) converges 
as z ^ 1. Consequently, we have Eq. (50). 
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